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The nonlinear amplitude modulation of known electrostatic plasma modes is examined in a generic 
manner, by applying a collisionless fluid model. Both cold (zero-temperature) and warm fluid de- 
scriptions are discussed and the results are compared. The moderately nonlinear oscillation regime is 
investigated by applying a multiple scale technique. The calculation leads to a Nonlinear Schrdinger- 
type Equation (NLSE), which describes the evolution of the slowly varying wave amplitude in time 
and space. The NLSE admits localized envelope (solitary wave) solutions of bright- (pulses) or dark- 
(holes, voids) type, whose characteristics (maximum amplitude, width) depend on intrinsic plasma 
parameters. Effects like amplitude perturbation obliqueness, finite temperature and defect (dust) 
concetration are explicitly considered. The relevance with similar highly localized modulated wave 
structures observed during recent satellite missions is discussed. 
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I. INTRODUCTION 



The amplitude modulation (AM) of waves is a generic nonlinear mechanism, which is long known to dominate 
finite amplitude wave propagation in dispersive media. In a generic context, the occurrence of AM is manifested as a 
slow variation of the wave's amplitude in space and time, which may be due to parametric wave coupling, interaction 
between high- and low- frequency modes or self- interaction of the carrier wave (auto- or seZ/-modulation) . The relation 
of this phenomenon to effects like secondary harmonic generation and modulational instability (MI) , possibly resulting 
in energy localization via localized pulse formation, is now long established in fields as diverse as Condensed Matter 
Physics, Nonlinear Optics and Biophysics Ql 0, 0] • With respect to plasma modes 5, :6], the occurrence of AM and 
MI has been confirmed by experiments related to the nonlinear propagation of electrostatic (ES, e.g. ion-acoustic) 
0, El as we U as electromagnet ic ( EM, e.g. whistler) waves. Early numerical simulations of electron cyclotron 
waves also predict such a behaviour 11]. 
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FIG. 1: Localized envelope structures in the magnetosphere (reprinted from [T^ l. 
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FIG. 2: Modulated structures, related to 'chorus' (EM) emission in the magnetosphere (CLUSTER satellite data; reprinted 
from H). 

In the context of Space Physics, localized modulated wave packets are encountered in abundance e.g. in the Earth's 
magnetosphere, in fact associated with localized field and density variations which were observed during recent satellite 
missions dllll 13 ; see e.g. Figs. H [JJ The occurrence of such wave forms is, for instance, thought to be related 
to the broadband electrostatic noise (BEN) encountered in the auroral region [lj- Furthermore, recent studies have 
supplied evidence for the relevance of such effects in dust-contaminated plasmas [Dusty or Complex Plasmas), where a 
strong presence of mesoscopic, massive, charged dust grains strongly affects the nonlinear and dispersive characteristics 
plasma as a medium |l5l Il6j . The modification of the plasma response due to the presence of the dust gives rise to 
new ES/EM modes, whose self-modulation was recently shown to lead to MI and soliton formation; these include 
e.g. the dust-acoustic (DA) [HE, [SHI and dust-ion acoustic (DIAW) ES modes d IM |H , in addition to 
magnetized plasma modes, e.g. the Rao EM dust mode |24j . 
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FIG. 3: Electrostatic noise wave forms, related to modulated electron-acoustic waves (FAST satellite data; figure reprinted 
from [121 )• The co-existence of a high (carrier) and a low (modulated envelope) frequencies is clearly reflected in the Fourier 
spectrum, in the right. 

The purpose of this brief study is to suggest a generic methodological framework for the study of the nonlinear 
(self-) modulation of the amplitude of electrostatic (ES) plasma modes. The general results provided in the following 
are valid for various ES modes. The generic character of the nonlinear behaviour of these modes is emphasized, so 
focusin g up on a specific mode is avoided on purpose. Where appropriate, details may be sought in the references 
[T^.l22ll23Ll25l l2r^. where some of this material was first presented. 

II. THE MODEL: FORMULATION AND ANALYSIS 

In a general manner, several known electrostatic plasma modes 0, are modeled as propagating oscillations 
related to one dynamical plasma constituent, say a (mass m a , charge q a = s a Z a e; e is the absolute electron charge; 
s = s a = (fa /| 9a | — ±1 is the charge sign), against a background of one (or more) constituent (s), say a' (mass 
?7ia', charge q a i = s a iZ a >e, similarly). The background species is (are) often assumed to obey a known distribution, 
e.g. to be in a fixed (uniform) or in a thermalized (Maxwellian) state, for simplicity, depending on the particular 
aspects (e.g. frequency scales) of the physical system considered. For instance, the ion-acoustic (IA) mode refers to 
ions (a = i) oscillating against a Maxwellian electron background (a' = e) 0,123, the electron- acoustic (EA) mode 
[H l2q] refers to electron oscillations (a — e) against a fixed ion background (a' = i), and so forth 0,0. As regards 
dusty plasma modes, the DA mode describes oscillations of dust grains (a — d) against a Maxwellian electron and 
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ion background (a' = e, i) while DIA waves denote IA oscillations in the presence of inertial dust in the 

background (a — i, a' = e, d) |l6l 122. l23|. Finally, this formalism readily applies in the case when a co-existence 
of two different populations of the same particle species occurs in the background, e.g. when two different electron 
temperatures are present (a' — c, h, for cold and hot electrons), affecting IA oscillations (a = i) |25| : this situation is 
witnessed in the upper parts of the Earth's atmosphere. 



A. A generic fluid description 

A standard (single) fluid model for the dynamic species a consists of the moment evolution equations: 

-7^- + V • (n a UoJ = 

— - + u Q -Vu a = -— V$ V Pa 

ot m a m a n a 

— ^+u a -Vp a = -7p Q V-u Q , (1) 
ot 

where n a , u a and p a respectively denote the density, mean (fluid) velocity and pressure of species a. The electric 
potential $ obeys Poisson's eq.: 

V 2 $ = -47r V" q a "n a » = Aire (n e - Z { n { + ...) (2) 

a"— a, {a'} 

Overall charge neutrality is assumed at equilibrium, i.e. q a n a Q = — Yl a ' 1 a ' n a'-0- The parameter 7 = cp/cy = 
1 + 2/ f denotes the specific heat ratio (for / degrees of freedom). 

By choosing appropriate scales for all quantities, the above system may be reduced to the following form: 

l + v-H = 0, (3) 

-^+u-Vu = -sV0-~Vp, (4) 
at n 

dp 

— +u-Vp = -7pV-u (5) 
ot 

(the index a will be understood where omitted, viz. s = s a ). The re-scaled (dimensionless) dynamic variables are 
now: n = n a /n a fl, u = u Q /c», p — p a /n aj okBT a , an d <f> — \q a |$/(fcsT»), where n a ,o is the equilibrium density 
and c* = (fcsT^/mc) 1 / 2 is a characteristic (e.g. sound) velocity. Time and space are scaled over appropriately 
chosen scales to [e.g. Lo~ x a = (47rn a) o9a/ m a)~ 1/ ' 2 ] an< ^ r o = c *^o! T a is the fluid temperature, and T* is an effective 
temperature (related to the background considered), to be determined for each problem under consideration (fc^ is 
Boltzmann's constant). The temperature ratio T a /T„ is denoted by <j, in this warm model 27] (the so-called cold 
model is recovered for a = 0; see that Eq. (J5J then becomes obsolete). The Lorentz force term was omitted, since 
wave propagation along the external magnetic field is considered. The system is closed by Poisson's equation, which 
may now be expressed as 

V 2 = -s [n + ^vv/K,o?£«)] =-s(n-n). (6) 

a' 

Note that the neutralizing background (reduced) density 

n = - V* Ua qa = — y^ s a ,Z a <n a/ (7) 

, n a ,o q a Sq Z a n a fi e —f 

is a priori [40l | a function |4l| of the potential 0; furthermore, it depends on the physical parameters (e.g. background 
temperature, plasma density, defect concentration, ...) involved in a given problem. The calculation in the specific 
paradigm of IA waves is explicitly provided below, for clarity. 
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B. Weakly nonlinear oscillation regime 

What follows is essentially an implementation of the long known reductive perturbation technique |28l l29l I30L l3l| , 
which was first applied in the study of electron plasma 28] and electron-cyclotron 11] waves, more than three decades 
ago. 

Equations 101-01 and (JSJ form a system of evolution equations for the state vector 42] S = {n,u,p,<f>} which 
accepts a harmonic (electrostatic) wave solution in the form S = S exp[i(fcr — Lut)] + c.c. Once the amplitude of 
this wave becomes non-negligible, a nonlinear harmonic generation mechanism enters into play: this is the first 
signature of nonlinearity, which manifests its presence once a slight departure from the weak-amplitude (linear) 
domain occurs. In order to study the (amplitude) modulational stability profile of these electrostatic waves, we 
consider small deviations from the equilibrium state S^ ^ = (1,0, 1,0) T , viz. S = S^ ^ + eS^ + e 2 S( 2 ) + where 

e < 1 is a smallness parameter. We have assumed that 43] S^ n) = Ydl-cc S jf( X > T ) e a(kr ~^ (for j = 1, 2, d + 3 
^ the condition Sj™^ = Sjf holds, for reality). The wave amplitude is thus allowed to depend on the stretched 

(slow) coordinates X — e(x — A) t and T — e 2 1; the real variable A = duj(k)/dk x = v g denotes the wave's 
group velocity along the modulation direction |44[ x. The amplitude modulation direction is assumed oblique with 
respect to the (arbitrary) propagation direction|4q. expressed by the wave vector k = (k x , k y ) — (k cos 6, k sin 9). 
Accordingly, we set: d/dt — > d/dt — ev g d/dX + e 2 d/dT, d/dx — > d/dx + ed/dX (while d/dy remains unchanged) 
and V 2 — > V 2 + 2ed 2 /dxdX + e 2 d 2 /dX 2 , in all the above evolution equations. 
By expanding near (4k0, Poisson's eq. may formally be cast in the form 

V 2 = (j>-a4> 2 + a'<f> 3 -sj3(n-l), (8) 

where the exact form of the real coefficients a, a' and (3 (to be distinguished from the species indices above, obvi- 
ously) are to be determined exactly for any specific problem, and contain all the essential dependence on the plasma 
parameters. Note that the right-hand side in Eq. (JSJ cancels at equilibrium. 

a. A case study: ion-acoustic waves In order to make our method and notation clear, let us explicitly consider 
the simple case of ions (a — i and q a = qi = +Zie, i.e. s a — +1) oscillating against thermalized electrons (a' = e and 
la 1 = Qe = — e, i-e. s a i = — 1; n e = n e fle e ®^ kBT °^). Adopting the scaling denned above, and using the equilibrium 
neutrality condition n ei o = Zin e> o, it is a trivial exercise to cast Poisson's Eq. @ into the (reduced) form: 

V 2 = r /c») 2 [n - = -(n - , (9) 

where we took: to = r$/c* = lu~\ and ^ = T»/ {Z{T e ). Now, expanding near cf> w 0, viz. e^ w 1 + £4> + ^ 2 (f> 2 /2 + 
^ 3 4> 3 /6 + we have: 

V 2 (/.«^ + f¥/2 + ^ 3 /6-(«-l). (10) 

Finally, setting the temperature scale T* equal to T* = Z{T e , for convenience (viz. £ = 1^)|46|. one recovers exactly Eq. 
(|HJ) with a — —1/2, a' = 1/6 and (3 = 1. It may be noted that this case has been studied, both for parallel modulation 
(9 = 0), in Ref. [23, and for oblique modulation, in Refs. [3^; those results are recovered from the formulae below. 

b. Amplitude evolution equations By substituting into Eqs. @ - © & n( i © an d isolating distinct orders in e, 
one obtains a set of reduced evolution equations in the new variables. One is then left with the task of isolating 
orders in e" (i.e. n = 1, 2, ...) and successively solving for the harmonic amplitudes Sj\ . The calculation, particularly 
lengthy yet perfectly straightforward, can be found e.g. in 22] for IA (s = +1) and in [2(| for EA waves (s = —1). 

The first harmonic amplitudes are determined (to order ~ e 1 ) as 

-si±f^^- «« = «a = i P « (id 

p uj cos 9 ' tJsmO y 7 

in terms e.g. of the potential correction = ip, along with the dispersion relation to 2 = (3k 2 /(k 2 + 1) + jerk 2 . 
Furthermore, the amplitudes of the 2nd and 0th (constant) harmonic corrections are obtained in ~ e 2 ; the lengthy 
expressions are omitted here, for brevity |4*?| . 



C. The envelope (nonlinear Schrodinger) evolution equation 



The potential correction ip is found to obey a compatibility condition in the form of a nonlinear Schrddinger-type 
equation (NLSE) 

,i + P0 + CW« = O. (12, 
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Both the dispersion coefficient P, in fact related to the curvature of the dispersion curve as P = d 2 Lu/2dk 2 — 
[u/'(fc) cos 2 9 + u>'(k) sin 2 0/k\ /2, and the nonlinearity coefficient Q 1 which is due to carrier wave self-interaction, 
are functions of k, 9 and /3, as expected (in addition to a, a', for Q). The exact general expressions obtained j^H 
may be tailor fit to any given electrostatic plasma wave problem (via the form of the parameters a, a', (3), in view 
of a numerical investigation of the wave's amplitude dynamics (e.g. stability profile, wave localization; sec in the 
following) . 



III. AMPLITUDE STABILITY PROFILE 

It is known (see e.g. in Refs. 0,0, that the evolution of a wave whose amplitude obeys Eq. I|12|) depends on 
the coefficient product PQ, which may be investigated in terms of the physical parameters involved. To see this, first 
check that Eq. (|12[) supports the plane (Stokes') wave solution ip = ipo exp(iQ\ipo\ 2 T); the standard linear analysis 
consists in perturbing the amplitude by setting: ip = tpo + ei/'i.o cos (kX — ujT). One thus obtains the (perturbation) 
dispersion relation: Co 2 — P k 2 (Pk 2 — 2Q|-0i i o| 2 ) . One immediately sees that if PQ > 0, the amplitude ip is unstable 
for k < y/2Q/P\i/}± o\. If PQ < 0, the amplitude ip is stable to external perturbations. 

This formalism allows for a numerical investigation of the stability profile in terms of parameters e.g. like wavenum- 
ber k, (oblique) perturbation angle a, temperature T a , background plasma parameters etc. In figure we have 
depicted the region PQ < (PQ > 0) in black (white) color, for IA waves; see in Ref. [23| for details. 
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FIG. 4: The region of positive (negative) values of the product PQ are depicted in white (black), in the wavenumber k - 
modulation angle a plane. The first two plots refer to IA waves: a — (cold model); a = 0.05 (warm model). Similar for the 
latter two, but for DIA waves (see in the text); we have taken a negative dust density: /i = n e ,o / '(ZiUifi) = 0.5 (from Ref. |23]'). 
The dust presence strongly modifies the stability profile (rather enhancing instability here). 



IV. ENVELOPE EXCITATIONS 



The modulated (electrostatic potential) wave resulting from the above analysis is of the form 



ctpo cos(kr 



tut + O) + Q(e 2 ), where the slowly varying amplitude|48( ^o(X, T) and phase correction 0(X, T) (both real functions 
of {X, T}; see 36] for details) are determined by (solving) Eq. I|12|) for ip = i[)q exp(iO). The different types of solution 
thus obtained are summarized in the following. 

c. Bright-type envelope solitons For positive PQ, the carrier wave is modulationally unstable; it may either 
collapse, due to external perturbations, or lead to the formation of bright envelope modulated wavepackets, i.e. 
localized envelope pulses confining the carrier (see Fig. [SJl: 
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(13) 



36, 37], where v e is the envelope velocity; L and represent the pulse's spatial width and oscillation frequency (at 
rest), respectively. We note that L and ipo satisfy i^o = (2-P/Q) 1 / 2 = constant [in contrast with Korteweg-deVries 
(KdV) solitons, where L 2 ip Q — const, instead]. Also, the amplitude ipo is independent from the velocity v e here. 

d. Dark-type envelope solitons For PQ < 0, the carrier wave is modulationally stable and may propagate as a 
dark/grey envelope wavepackets, i.e. a propagating localized hole (a void) amidst a uniform wave energy region. The 
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FIG. 5: Bright type modulated wavepackets (for PQ > 0), for two different (arbitrary) sets of parameter values. 
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FIG. 6: Bright envelope soliton propagation, at different times ti < • ■ • < fcj (arbitrary parameter values): cf. the structures 
encountered in satellite observations, e.g. see Fig. 



exact expression for dark envelopes reads |3fJ, |3 



tanh 



[ X-v e T 
V U 



6= ±[ Ve X + (2PQ^ Q 2 ^f 



T 



(14) 



(see Fig. |7K); again, L'tp' = {2\P/Q\) 1 / 2 (=cst.). 

e. Grey-type envelope solitons The grey-type envelope (also obtained for PQ < 0) reads |36ll31 



V>o = ^" {l-d 2 scch 2 {[X-v e T]/L"}}^ 2 , 



G =2P 



iii-V ( -V 2 -2PQri)T + Q i 



— S sin 1 



d tanh 



1 — d 2 sech 2 



2 / X- 



-,1/2 



(15) 



Here 0o is a constant phase; S* denotes the product S = sign(P) x sign(w e — Vb). The pulse width L" = 
(\P/Q\) 1 / 2 /(difj" Q ) now also depends on the real parameter d, given by: d 2 = 1 + (v e — V ) 2 /(2PQi(>"l) < 1. 

Vq = const, e satisfies (21113 : ^o- \]2\PQ\i^l < v e < V + \j 1 2\PQ\ijj"l. For d = 1 (thus V = u e ), one recovers 
the dark envelope soliton. 





FIG. 7: Dark (left) and grey (right) type modulated wavepacket (for PQ < 0). See that the amplitude never reaches zero in 
the latter case. 



V. CONCLUSION 



This work has been devoted to the study of the conditions for occurrence of modulational instability, related to the 
formation of envelope localized structures, with respect to electrostatic waves propagating in an unmagnetized plasma. 
We have shown that the envelope modulated electrostatic wave packets which are widely observed during satellite 
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missions and laboratory experiments, may be efficiently modeled by making use of a reductive perturbation (multiple 
scales) technique |28j . Explicit criteria are thus obtained, which determine the wave's modulational stability profile 
and predict the occurrence of localized envelope excitations of either bright or dark/grey type. This methodology 
allows for an investigation of the nonlinear modulational profile of a (any) given electrostatic mode. Generalization 
in the presence of a magnetic field is on the way and will be reported soon. 
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